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ABSTRACT 

The  purpose  of  the  InYestlg&tion  was  to  determine  the 
feasibility  of  obtaining  an  electrical  circuit  whose 
behavior  would  be  a  satisfactory  analog  of  the  dynamics 
of  a  missile  in  flight.  It  was  originally  hoped  that  this 
could  be  done  with  passive  circuit  elements;  however,  it 
was  found  that  due  to  severe  non-linearity  and  unilateral 
coupling,  it  was  not  possible  to  establish  a  passive  circuit 
to  represent  even  a  linear  approximation  of  the  missile 
system. 

It  was  also  hoped  to  find  a  mathematical  procedure  for 
analyzing  the  analog  circuit  which  would  allow  analysis  of 
the  missile  system  by  analogous  responses.  Due  to  the 
severity  of  the  non-linearity  this  was  not  possible.  Because 
of  the  nature  and  Importance  of  the  stability  problem  in 
non-linear  systems,  it  was  recommended  that  the  analysis  of 
the  missile  system  be  accomplished  by  utilizing  an  active 
electrical  analog  simulator  arranged  specifically  for  solution 
of  the  set  of  behavior  equations  representing  missile  dynamics 
including  wind  profile  forcing  functions.  A  study  was  made 
of  the  various  devices  which  could  be  used  to  perform  the  cir¬ 
cuit  operations  Involved.  This  study  lead  to  the  belief  that 
the  most  practical  simulator  would  be  composed  of  standard 
analog  components. 
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SECTION  I 
INTRODUCTION 

This  report  covers  research  conducted  under  Contract 
DA-36-03U-0RD-3513RD  between  AOMC  and  The  Pennsylvania  State 
University,  The  work  was  done  in  the  Department  of  Electri¬ 
cal  Engineering. 

1.1  Object 

The  object  of  the  investigation  was  to  determine  to  what 
degree  it  would  be  possible  to  utilize  electric  circuits  as 
analog  devices  so  arranged  as  to  simulate  dynamic  missile 
behavior, 

1.2  Need  for  the  Study 

In  the  ideal  case  the  provision  of  an  electrical  analog  of 
the  missile  system  would  consist  of  representing  the  para¬ 
meters  of  the  missile  system  in  terms  of  electrical  para¬ 
meters,  This  would  permit  personnel  trained  primarily  in 
electrical  theory  to  analyze  the  missile  system  in  terms  of 
electrical  parameters  with  which  they  are  more  acutely 
familiar.  In  utilizing  the  electrical  analog  the  procedure 
would  be  to  obtain  a  set  of  behavior  equations  for  the 
electrical  system.  This  would  consist  of  a  set  of  differ¬ 
ential  equations  the  solution  of  which  would  be  the  response 
of  the  electrical  system.  This  response  could  then  be 
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interpreted  in  terns  of  the  parameters  and  response  Tariables 
of  the  missile  system  yielding  the  desired  dynamic  behaTior. 
The  success  of  this  procedure  depends  upon  obtaining  a  use- 
able  mathematical  solution  of  the  behavior  equations*  Any 
solution  which  is  not  in  closed  form  would  be  unsatisfactory 
because  of  the  obscurity  of  the  information  desired*  What 
is  needed  is  a  method  which  can  be  applied  to  the  system  of 
equations  in  a  manner  similar  to  the  way  in  which  the  Laplace 
transform  Is  applied  in  obtaining  the  solution  to  systems  of 
linear  equations* 

In  the  absence  of  a  satisfactory  mathematical  solution  an 
alternative  procedure  would  be  to  provide  the  physical 
electrical  analog  system  and  interpret  its  response  in  terms 
of  the  response  of  the  missile  system* 
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SECTION  II 
ANALYSIS 


2,1  General  Prooedtir> 

The  first  step  In  approaching  the  problem  was  to  obtain  from 
AOMC  a  set  of  equations  representing  the  dynamic  behavior  of 
a  missile  system.  These  equations  are  shown  on  Page  5  of 
reference  15,  Supplementary  to  this  the  limited  literature 
search,  references  2,  3,  4,  5,  6,  7,  12,  16,  was  conducted 
to  determine  whether  these  equations  were  presented  in  the 
best  form  and  to  gain  familiarity  with  the  general  approach  to 
the  study  of  missile  dynamics.  The  result  of  this  revealed 
that  the  equations  as  presented  in  reference  1$  are  sub¬ 
stantially  complete  and  representative  of  a  practical  missile 
system.  These  equations  were  used  ae  the  basis  for  analysis 
throughout  this  report. 

Having  obtained  a  satisfactory  set  of  equations,  these  were 
examined  to  establish  the  mechanical  parameters  involved  and 
the  appropriate  variables  normally  employed  in  missile 
dynamics  studies,  A  study  was  then  made  to  determine  an 
appropriate  set  of  electrical  parameters  and  variables  repre¬ 
senting  a  suitable  analog  system.  Since  no  analog  system  is 
necessarily  unique  then  many  analogs  may  exist,  it  was 
necessary  to  defer  using  a  specific  system  until  a  study  could 
be  made  of  the  mathematical  operations  and  hence  of  the  electric 
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circuit  operations  involved  in  representing  each  tern  of 
each  equation  in  the  set  of  behavior  equations. 

Having  established  circuit  operations  involved  in  the 
analog  system  a  study  was  then  made  of  the  various  electrical 
devices  and  circuits  capable  of  performing  these  operations. 

Examination  of  the  behavior  equations  taken  from  reference  1$ 
as  shown  in  section  2.2  shows  that  they  are  severly  non-linear. 
In  view  of  this  it  wac  decided  to  obtain  a  first  approximation 
to  these  equations  by  applying  a  linearizing  technique.  An 
analog  circuit  was  then  developed  from  the  linearized  equations 
and  was  studied  to  determine  the  best  choice  of  circuit  vari¬ 
ables  and  parameters.  It  was  immediately  discovered  that  the 
original  equations  exhibit  the  characteristic  of  unilateral 
coupling.  This  implies  that  the  electrical  analog  circuit 
must  therefore  contain  unilateral  devices.  It  was  also  found 
that  •isolation”  is  a  circuit  operation  which  must  be  provided. 

After  the  study  of  the  linear  circuit  was  completed  a  new 
analog  circuit  was  developed  embodying  all  of  the  non-linear 
terms  of  the  original  set  of  equations.  This  circuit  was 
studied  and  conclusions  drawn,  details  of  which  are  presented 
later. 
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circuit  operations  Involved  In  representing  each  term  of 
each  equation  In  the  set  of  behavior  equations* 

Having  established  circuit  operations  Involved  In  the 
analog  system  a  study  was  then  made  of  the  various  electrical 
devices  and  circuits  capable  of  performing  these  operations* 
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2,2  The  Equatlona 
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Consider  first  the  set  of  equations  taken  from  reference  1$ 
which  are  based  on  Figure  1  taken  from  the  same  reference. 


The  nominclature  used  in  the  discussion  of  these  equations 
is  as  follows: 
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the  moment  of  Inertia 

the  coefficient  of  lift 

the  diameter  of  the  rocket 

the  density  of  the  air 

the  geometrical  yeloeity 

the  path  angle  (V  local  rertical) 

the  attitude  angle  (long,  axis  ^  local  yertlcal) 

gyro  program  angle  (vertical  at  launch) 

the  angle  of  attack  (axis  flow  direction) 

the  wind  angle  (V  V^,) 

control  deflection  angle 

the  angle  at  earth  center  between  launch  and 
present  position 

the  mass  of  the  rocket 

gravitation 

the  thrust  in  direction  of  axis 
the  axial  air  force 
the  normal  air  force 
the  normal  vane  force 
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V  ■  the  Telocity  relative  to  surrounding  air 
r 

a  -  ■  the  coefficient  of  linear  terms  on  B 
0 

h  m  the  coefficient  of  a  terms  In  6 
0 

■  the  restoring  moMnt  per  unit  moment  of  Inertia 
C2  ■  the  control  moment  per  unit  moment  of  inertia 

d  ■  the  damping  moment  coefficient  per  unit  moment  of 

Inertia > 

The  forces  in  the  direction  of  the  missile  Telocity,  divided 
by  the  mass  of  the  missile  give  the  equation 

V  •  cos(a-ay)  -  sin(o-o^)  -  g  cos  •©,  (2-1) 

while  the  forces  perpendicular  to  the  velocity,  divided  by 
MV  give  the  equation 

0  +  ^  ■  8ln(o-oi^)  +  co8(a-o^)  +  ^  sin  6.  (2-2) 

The  moments  about  the  center  of  gravity  of  the  missile 
divided  by  the  moment  of  inertia  yield  the  equation 

(/.+  d  ♦  C^a  +  CgP  •  0 
The  idealized  control  equation  becomes 

P  •  *0^^*  l^-x)  +  »x(C^+^^  -x)  + 
and  from  angular  relations 

•  “  ^  ^2-5) 


(2-3) 


(2-4) 
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It  ia  realized  that  many  of  the  parameters  are  not  constant 
and  therefore  constitute  either  time  dependent  coefficients 
or  non-linear  terms,  but  for  the  present  work  It  was  decided 
to  consider  only  the  following  fire  rariables. 


V  ■  the  velocity 
a  ■  the  attack  angle 
|EI  ■  the  center  angle 
0  ■  the  path  angle 
(j^»  the  attitude  angle 


With  the  use  of  the  relations  found  In  reference  15,  any  one 
of  these  five  variables  can  be  removed  by  direct  substitution* 
This  leaves  four  variables  and  only  three  equations*  A  fourth 
equation  can  be  derived  using  trigometric  relationships  from 
Fig.  1. 


Since  •  V  sin  0,  V  •V  cos  0,  and  y 

X  y  r  • 

then  y  ■  cos  0  dt. 

Now,  X  ■  (Rq  +  y)  0 

from  which  (R^  +  y)  0  ■  V  sin  0 

which,  when  differentiated  yields 


V  sin  0  •  (R^  ♦  y)  0  +  y  0» 
From  which 


V  sin  0  ■  [«,  */-  sin  0  d^  ^  +  V(cos  C)  0  (2-6) 
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utilizing  this  fourth  equation,  we  now  have  a  set  of  four 
sinultaneous  integro-differential  equations  describing,  to 
some  degree  of  approximation,  the  dynamic  response  of  a 
missile.  If  these  equations  were  linear,  well  known  tech¬ 
niques  for  handling  linear  sets  could  be  applied  and  the 
entire  missile  behavior  could  readily  be  obtained. 

Examination  of  the  equations  will  reveal  that  they  are 
decldely  non-linear. 

It  is  possible  to  solve  these  equations  by  numerical  methods. 
However,  such  a  procedure  would  best  be  done  by  employing  a 
digital  computer.  Although  this  would  yield  a  solution,  this 
procedure  would  be  contrary  to  the  objective  of  this  research. 
It  was  felt  desirable  to  consider  as  a  first  approximation  a 
linear  model  obtained  by  applying  a  linearizing  technique  to 
the  equations. 

2.3  The  Linear  Model 

As  a  first'  approximation  to  the  physical  system,  it  was  de¬ 
cided  to  linearize  equations  (2-1)  through  (2-5).  This  was 
done  in  order  to  provide  the  basis  for  a  linear  electrical 
network  whose  behavior  would  approximate  the  behavior  of  the 
missile. 

Using  the  first  term  of  the  Maclaurin  series  for  sine  and 
cosine  functions,  equation  (2-1)  becomes 
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V  •  -  +  -Sjjr  •  it 

and  equation  (2-2)  becomea 

h.i 

Substituting  the  angular  relation  (2-5)  in  these  equations 
yields 


and 


;  N+R  //.  P-I 

V  -  -  (  Cf-  0)  +  -TT  “  «* 


0  ♦  ^  -  Iji  (6f-  «)  +  +  f  ft. 


(2-7) 


(2-8) 


Substituting  equations  (2-U)  and  (2-5)  into  equation  (2-3) 
and  ignoring  the  higher  derlTative  terms  in  (2-U))  yields 


For  convenience,  equations  (2-6),  (2-7)  and  (2-8)  can  bo 
rewritten 


(2-9) 


7  •  a 


i<(?- «) 


•  ♦  Ag 


0  ♦  ^  ^  (ce-  0)  -  ^  ‘ 


1-0 


ie  *  ♦  *6^^“  «)  •  A3  +  A^^ 


(2-10) 

(2-11) 

(2-12) 


where 


‘1  " 


NfR 
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*2 

*3 

•u 

»6 


g 
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^2*^0 


•re  considered  constant,  and 


Ag  ■  g 

^3  *  *0^2^ 


can  be  considered  forcing  functions  on  the  system*  Eqnatlons 
(2-10)  through  (2-12)  are  still  non-llnear  because  of  the 
existence  of  terms  containing  two  warlables*  They  can  best 
be  made  linear  by  considering  the  mlselle  velocity,  V,  to  be 
a  constant.  This  will  yield  results  valid  only  over  a  short 
excursion  of  flight,  iifhlle  this  Is  a  serious  restriction  on 
the  accuracy  of  the  problem,  it  must  be  remembered  that  this 
linear  set  of  equations  is  intended  only  as  a  first  approxi¬ 
mation. 


With  the  velocity  a  constant,  equations  (2-9)  through  (2-11) 
become 

«l(C^-  0)  -  Aj  ♦  Ag  (2-13) 
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(2-llt) 


Q  +  0  -  0)  -  aj  e  •  Aj 

^  ^5* 

These  are  linear  differential  equations  in  ^ ,  0,  and 
where : 


(2-15) 


a,  • 


N+R 


*2 

*3 

\ 

*5 

»6 


F-I 

w 

I 
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°2*o 


^1  *  °2^o» 


are  considered  constants,  and 


4  F-I 

^1  ■  “IT 


Ag  ■  g 
»  N+R 

\  -  »o°2X 

A5  ■  -(Cl  +  ^o°2^S 


can  be  considered  forcing  functions  on  the  system. 


Based  on  the  linearized  equations  (2-13)  through  (2-15),  an 
effort  was  made  to  obtain  a  linear  passlre  electrical  analog 
circuit,  having  behavior  equations  of  the  same  form.  This 
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lamediately  to  difficulty  because  the  coupling  terns 
rereal  a  non-bilateral  charncterlstlc  uhlch  Is  not  attain¬ 
able  by  passive  linear  circuit  elenents.  This  leads  to  the 
conclusion  that  there  is  no  passive  linear  multl-mesh 
electrical  network,  whose  behavior  can  be  expressed  by- 
equations  (2-13),  (2-ll()  and  (2-l^)»  It  should  be  noted 
that  although  these  equations  are  linear,  the  coefficients 
of  the  coupling  terns  are  not  equal. 

The  unilateral  coupling  exhibited  In  these  equations  requires 
the  introduction  of  unilateral  devices  Into  the  analogous 
electrical  network.  The  characteristics  of  such  devices  nust 
be  carefully  defined.  The  unilateral  device  is  synbolically 
shown  In  Figure  2,  lere  the  device  is  represented  as  a 
rectangle  having  two  input  terninals  and  two  output  terminals. 
The  Ideal  device  is  assumed  to  have  infinite  input  Impedance 
and  zero  output  Impedance  and  possesses  no  common  ground.  The 
capital  letter  in  the  upper  center  of  the  rectangle  indicates 
the  function  being  performed.  For  the  linear  model  there  are 
only  three  such  functions,  namely  inter gratlon  represented  by 
a  capital  I,  differentiation  represented  by  a  capital  D,  and 
amplification  represented  by  a  capital  A,  The  small  letter  in 
the  lower  center  of  the  rectangle  represents  the  scaler 
multiplier  associated  with  the  function.  Since  the  device  is 
a  four  terminal  device  with  input  and  output  terminal  pairs 
Isolated,  it  is  possible  for  the  output  to  be  connected  in 
either  of  two  possible  arrangements  to  provide  the  desired 
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polarity*  The  dot  aeeoclated  irith  Inpvt  and  ontptit  terminals 
indicates  'the  instantaneous  polarity*  That  is,  when  the  in¬ 
put  dotted  terminal  is  driven  positive  the  output  dotted 
terminal  will  be  positive* 

It  was  initially  decided  that  in  the  linear  electrical  analog 
circuit  the  various  mechanical  quantities  should  be  represent¬ 
ed  by  currents  in  which  case  the  unilateral  coupling  will  re¬ 
sult  in  a  potential  introduced  in  one  loop  by  vlrture  of  a 
current  in  a  different  loop.  With  this  in  mind.  Figure  3  shows 
how  the  device  can  be  used  to  provide  unilateral  coupling* 
Cqpsider  that  is  d  loop  current  in  loop  1  which  passes 
through  resistor  and  it  is  desired  to  Induce  a  potential 
in  loop  No*  2  which  is  proportional  to  this  current  while  at 
the  same  time  prohibiting  I2  in  loop  2  from  inducing  any  po¬ 
tential  in  loop  1*  The  input  terminals  of  the  unilateral  de¬ 
vice  are  connected  across  Rj^*  Since  the  input  Impedance  of 
the  device  is  infinite,  the  device  offers  no  loading  on  loop 
No.  1.  The  output  voltage  induced  in  loop  2  will  then  be 
described  by  the  equation  in  the  figure.  It  will  be  noted  that 
the  polarity  marks  indicate  that  the  potential  induced  in  loop 
2  is  positive  for  a  positive  current  in  loop  1*  Had  it  been 
desired  that  this  coupling  be  negative,  it  would  only  have  been 
necess^ary  to  reverse  the  output  terminals  or  the  input  ter¬ 
minals  but  not  both. 

The  resistance  R^  as  shown  in  Figure  3  not  only  provides  the 


Fig.2 
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potential  for  unilateral  ooupling  to  loop  2  but  also  provides 
a  self-lmpedanoe  for  loop  1*  Since  the  function  provided. by 
the  device  Is  not  always  anpllf loatlon  but  may  be  Intergratlon 
or  differentiation,  Figures  li  and  5  show  how  the  coupling 
could  be  arranged  to  provide  potential  In  loop  2  which  Is  the 
Integral  of  the  current  In  loop  1.  It  will  be  noted  that 
this  operation  Is  not  unique  but  that  such  an  operation  could 
be  provided  by  either  the  circuit  of  Figure  It  or  that  of 
Figure  5*  The  difference  between  these  two  circuits  however, 
lies  In  the  self-impedance  component  which  Is  desired  In  loop 
!•  In  Figure  U  the  self -Impedance  Is  resistance  and  the  In¬ 
tegration  Is  performed  by  the  device.  In  Figure  5  the  self- 
Impedance  is  that  of  the  capacitor  which  Integrates  within 
Itself,  thus  the  device  merely  amplifies  to  provide  the  desired 
operation.  It  is  evident  that  differentiation  could  be 
accomplished  by  using  an  inductor  In  loop  1  and  permitting  the 
device  to  simply  amplify.  On  the  other  hand,  if  the  self¬ 
impedance  desired  Is  resistance,  then  the  device  would  be  re¬ 
quired  to  differentiate. 

Using  these  basic  operations^  consider  the  circuit  of  Figure  6 
and  the  three  linear  equations  (2-13)  through  (2-l5)  which  are 
presented  by  the  three  Independent  circuit  loops.  Loop  No.  1 
represents  equation  (2-13),  loop  No.  2  equation  (2-ll)),  and 
loop  3  equation  (2-15).  Let  the  variable  in  loop  1  be  0  and 
be  represented  by  I^.  Then  the  first  term  of  equation  (2-13) 

Is  represented  by  a  potential  drop  across  the  resistance  which 
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l8  the  total  resistance  in  the  loop.  The  drive  functions 
and  A2  are  shown  as  generators,  with  polarities  as 
indicated.  The  second  tern  of  equation  (2-13)  is  a  potential 
equal  to  a  constant  multiplier  A^  times  the  variable  te  which 
is  the  integral  of  current  in  loop  3  representing  • 

This  potential  is  provided  by  unilateral  device  No.  1  shown 
at  the  left  where  the  input  potential  taken  across  the 
capacitor  in  loop  3  is  multiplied  by  the  constant  a^  and  is 
induced  into  loop  1  in  the  negative  direction.  Since  the 
variable  in  loop  3  is  ^  is  obtained  by  the  Integrating 

properties  of  the  capacitor.  It  will  be  observed  that  there 
are  no  other  unilateral  devices  feeding  into  loop  1  although 
there  are  three  others  feeding  from  loop  1  to  other  loops. 

In  equation  (2-lU),  represented  by  loop  2,  let  it  be  observed 
that  I2  represents  the  variable  0  even  though  0  Itself  does 
not  appear  in  the  equation.  This  choice  was  made  for  conven¬ 
ience  in  obtaining  the  necessary  coupling  terms  to  and  from 
other  loops.  The  first  term  of  equation  (2-11;),  0  ,  is  the 
potential  across  the  total  Inductance  in  loop  2.  The  second 
term  0  is  obtained  by  unilateral  device  No.  5  which 
differentiates  the  potential  developed  across  a  resistance 
in  loop  1.  The  value  of  the  resistance  in  loop  1  is  not 
critical  as  long  as  the  product  of  this  value  of  resistance 
and  the  gain  function  of  the  unilateral  device  is  equal  to 
unity.  Generally  in  such  cases  the  resistance  would  be  chosen 
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LOOP  3 


FIGURE  6 


as  unltx*  Tha  third  tera  of  aquation  (2-lU)  is  a  constant 
nultlpllar  tlaas  Q  which  is  obtalnod  by  unilateral  derioa 
No*  6*  Assuaing  tha  rasistanea  in  loop  1  is  chosen  as  unity 
than  tha  gain  function  of  tha  unilateral  dsTlca  is  set  at 
tha  number  (a2  -  a^)*  It  will  be  noted  that  the  polarity  is 
arranged  to  provide  a  positive  potential  in  loop  2  due  to  a 
positive  current  in  loop  1* 

The  fourth  term  of  equation  (2-lU)  is  a  constant  multiplier 
times  the  variable  and  should  appear  in  the  negative 
direction*  This  is  obtained  from  unilateral  device  No*  3 
which  simply  amplifies  the  potential  taken  across  a  capacitor 
in  loop  3  producing  a  voltage  in  loop  2  equal  to  the  number 
a2  multiplied  by  the  integral  of  which  is  6^*  The  drive 
function  A^  is  again  shown  as  a  generator*  Equation  (2-l5), 
represented  by  loop  3t  can  be  analyzed  in  the  sane  manner  where 
it  will  be  noted  unilateral  device  No*  2  provides  the  fifth 
term  and  device  No*  U  provides  the  fourth  term* 

The  diagram  shown  in  Figure  6  is  not  unique  and  other  such 
diagrams  could  be  arranged  to  represent  the  same  set  of 
equations*  The  one  shown  has  been  devised  as  an  example  to 
show  the  application  of  the  unilateral  devices  and  the 
flexibility  in  choice  of  variables  represented  by  currents  in 
the  various  loops* 

The  behavior  equations  of  this  circuit  are  of  the  form  of 
equations  (2-13 }  through  (2-l5)«  and  in  terns  of  the  electrical 
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paraneters  are  as  follove: 


♦  Ag  (2-16) 

Lglg  +  +  (»2  “  ■  ^3  (2-17) 

li^l^  +  ^  ♦  a^lg  +  a^i^^  ■  A|^+A^  (2—18) 


where t 

R 

1 

L 

1 

L 

R 

C 

i 


1 

1 

2 

2 

3 

3 

3 

3 


*1  from  equations  (2-13)  and  (2-l6) 

0 


from  equations  (2-lU)  and  (2-1?) 


1 

*u 

l/(a^4ag)  from  equations  (2-l5)  and  (2-18) 


A 

This  shows  the  relationship  between  the  mechanical  and 
electrical  parameters* 


Since  the  circuit  of  Figure  6  represents  only  a  first 
approximation,  it  is  now  desirable  to  obtain  a  circuit 
which  more  accurately  represents  the  missile  system.  This 
is  accomplished  by  utilizing  the  non-linear  equations  without 
simplification. 


2.1i  The  Non-Linear  Model 

If  the  mass,  M  along  with  0,  0,  ^ ,  and  V,  is  allowed  to  vary. 
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equatlona  (2-1),  (2^-2),  (2-3)  and  (2-6)  can  be  wrlttenr 


V  >  ^  cos  (d^- 

0)  -  -1  sin  (^-  e)  -  cos  e 

(2-19) 

0  +  fl  ■  ^  ■!» 

(C^-  6)  +  ^  coe  (^-  0) 

4-  ^  sin  0 

(2-20) 

^  ♦  D  ♦  B  . 

.  FO4-O^4-H0>  Q(t) 

(2-21) 

T  eln  «  -  J  0 

/  V  cos  fi  dt  ♦  V(cos  0)  li 

(2-22) 

where: 

A  -  F-X 

6  -  N-t-R 
C  -  g 

D 
E 
F 
0 


d  ♦  ^2^1 

*  °2*o  *  ‘^2'’o 
^ 


^2*0 


H  -  Cga^^ 


R 


Q(t)  Drive  •  -(C^+Cgb^)  +  Cga^^X  +  Cga^^X 


(2-23) 


In  the  equation  of  the  forcing  function^Q(t )^  (C^+C2b^)o^ 
represents  the  wind,  and  Cga^X  +  Cga^^I  represents  the  pre¬ 
programmed  gyro  setting* 
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A  four  B«8h  eleotrlcal  analog  network  was  developed  without 
benefit  of  any  approximation  beyond  that  already  contained 
in  the  equations*  This  circuit  is  necessarily  complicated 
because  it  contains  all  of  the  non-linear  terms.  Figure  7 
should  be  compared  with  Figure  6* 

In  Figure  7  each  rectangle  represents  a  mathematical  circuit 
operation*  There  are  six  basic  operations  performed  by  the 
following  devices? 

M  ■  multiplier 
R  -  resolver 
A  ■  amplifier 
I  ■  integrator 
S  •  summer 

0  •  inverter  (reciprocal  generation) 

Each  multiplier,  designated  M,  has  an  output  equal  to  the 
product  of  its  two  Inputs.  Each  resolver,  designated  R,  has 
two  output  terminals  providing  simultaneous  sine  and  cosine 
functions  of  the  input  signal.  The  output  of  the  summer, 
designated  S,  is  the  sum  of  its  two  input  signals.  The  out¬ 
put  of  the  Inverter,  designated  Q,  is  the  inverse  of  its 
input.  Finally  the  amplifier  and  Integrator  are  shown  with 
double  letter  notation  where  the  upper  letter  represents  the 
function  (A  in  the  case  of  amplifier  and  I  in  the  case  of 
integrator)  and  the  lower  letter  represents  the  gain  which 
can  be  either  positive  or  negative.  Since  each  of  these 
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rectangles  represents  a  derice  performing  a  mathematical 
operation  it  is  assumed  that  a  suitable  device  can  be  provided 
which  will  perform  this  operation  to  a  satisfactory  degree* 

Such  devices  will  necessarily  have  to  be  unilateral  within 
themselves,  thus  the  arrowheads  on  the  diagram  show  the 
direction  of  signal  flow* 

Basically  the  circuit  of  Figure  7  is  a  four  mesh  electrical 
network*  These  meshes  are  shown  in  the  corners  of  the  dia¬ 
gram  and  are  numbered  1  through  1|*  They  correspond  to  the 
variables  in  equations  (2-19)  through  (2-22),  For  the  sake 
of  simplicity,  the  diagram  is  drawn  as  a  single  line  diagram 
except  at  the  points  where  signals  leave  or  enter  the  meshes. 
For  example,  consider  mesh  1  in  which  the  variable  is  velocity* 
The  only  signal  leaving  the  mesh  is  V  which  is  the  potential 
across  the  inductor  fed  to  an  integrator  over  two  leads*  The 
output  of  that  Integrator  is  a  single  lead  feeding  other  loops* 
It  is  clear  that  the  Integration  cancels  the  differentiation 
produced  by  the  inductor*  Hence,  the  output  of  the  integrator 
is  essentially  the  sum  of  the  voltages  entering  the  mesh* 
Although  this  could  have  been  obtained  by  a  simple  summer,  it 
was  desired  at  this  point  to  retain  the  identity  of  the  mesh 
representing  equation  (2-19)*  It  will  be  shown  in  section  3 
that  this  is  neither  necessary  nor  desirable* 

All  voltages  entering  the  mesh  are  shown  entering  in  a  series 
of  double  lead  circuits  frmm  the  amplifiers  shown  immediately 
to  the  right  of  the  mesh*  In  each  mesh  the  number  associated 
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with  the  eonponent  (that  is  inductor,  resiator,  or  capacitor) 
represents  the  coefficient  in  the  respective  equation.  It 
should  be  possible  now  to  follow  the  signals  through  the 
remainder  of  the  diagram  and  to  observe  the  generation  of 
various  non-linear  and  unilateral  coupling  terms. 

It  will  be  noted  that  the  currents  in  meshes  1  and  It  represent 
the  variables  V  and  0  respectlvelj,  while  in  meshes  2  and  3 
the  currents  represent  the  variables  0  and  respectively. 
The  choice  of  currents  representing  derivatives  in  meshes  2 
and  3  was  made  to  avoid  the  use  of  differentiators  because  of 
the  practical  problems  Involved  in  their  operation.  Mesh  3 
contains  the  drive  on  the  system  which  is  shown  in  equation 
(2-20).  This  is  a  potential  generator  representing  wind  angle 
and  preprogrammed  gyro  setting.  Each  of  the  four  variables 
can  be  read  out  at  the  appropriate  terminal  shown  on  the  dia¬ 
gram  as  a  small  circle.  Initial  conditions  are  introduced  on 
each  of  the  integrators  and  as  currents  in  loops  1,  3,  and  U. 
The  circuits  providing  these  initial  conditions  are  not  shown 
because  they  are  of  secondary  importance. 

A  careful  study  of  this  diagram  will  reveal  the  importance  of 
each  functional  block  and  the  part  it  plays  in  establishing 
the  necessary  intermesh  coupling.  It  is  expedient  at  this 
point  to  refer  to  Fig.  8  which  shows  a  basic  analog  simulator 
diagram  for  the  solution  of  the  same  h  simultaneous  equations. 
This  simulator  is  comprised  of  standard  analog  computer  com- 
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ANALOG  FUNCTIONAL  DIAGRAM  REPRESENTMG  THE  SYSTEM  OF 


pon«Dts  and  is  Introdneed  here  for  eonparison*  In  reriewlng 
Figure  7  the  reader  may  gain  the  impression  that  the  circuit 
is  unduly  complicated,  however,  comparison  with  Figure  8  will 
reveal  that  the  simulator  diagram  is  of  approximately  equal 
complexity.  In  Figure  7  there  are  32  functional  blocks,  in 
Figure  8  there  are  29  functional  blocks.  Considerable  time 
was  devoted  to  an  attempt  to  reduce  complexity  of  Figure  7 
but  with  the  support  of  the  circuit  of  Figure  8  it  was  con¬ 
cluded  that  very  little  simplification  is  possible.  It  is 
evident  that  most  of  the  diagram  is  comprised  of  functional 
blocks  which  are  necessitated  by  the  complex  non-linear  terms 
and  unilateral  coupling  terms.  This  appears  to  an  extent 
which  obscures  the  basic  four  mesh  concept.  Figure  7  was  pre¬ 
pared  utilizing  the  well  known  transfer  function  notation 
representing  the  behavior  of  the  individual  devices.  It  is 
assumed,  therefore,  that  each  device  performs  the  operation 
shown  Independent  of  the  load  on  its  output  terminals.  If 
this  is  not  the  case  and  loading  must  be  considered,  the 
determination  of  an  analog  network  becomes  many  times  more 
complicated,  if  indeed  it  is  at  all  possible. 

It  should  now  be  noted  that  although  the  circuit  of  Figure  7 
is  based  on  the  four  mesh  concept  where  currents  represent 
variables  in  the  equations,  all  the  functional  blocks  pro¬ 
viding  inter-mesh  coupling  are  potential  devices  which  have 
voltage  input  and  output.  It  is  not  necessary  or  unique  that 
the  variables  be  represented  by  currents  in  the  meshes.  It 
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is  possibl*  to  olinlnate  the  Beahaa  as  auch  by  replaeifig 
than  by  approprlata  additional  functional  bloeka.  For 
axanpla,  naah  No.  1  could  ba  aliniAatad  along  with  its  out¬ 
put  integrator  by  aimple  auauaing  of  the  voltagea  induced  into 
this  Kssh  and  letting  the  output  of  that  suBBer  be  the 
rarlabla  V.  SlBllarly  Bash  2  could  be  ellBinated  by  suBBing 
all  the  voltages  induced  in  that  Bssh  and  letting  the  output 
potential  ba  Q,  Meshes  3  and  U  could  alnilarly  be  ellninated^ 
howeveri  the  net  results  of  this  would  be  to  provide  a  cir¬ 
cuit  dlagran  oonprlsed  eonpletely  of  functional  blocks,  ob¬ 
scuring  the  concept  of  Bashes.  Such  a  dlagran  would  reduce 
to  that  of  Figure  8  where  all  functional  blocks  are  perforned 
by  standard  sinulator  elenents.  It  is  not  surprising  that 
this  should  happen  since  analysis  of  circuits  of  this  nature 
has  been  the  notlvating  influence  leading  to  the  developnent 
of  the  walog  conputer. 

2 .5  Stability  Considerations 

It  was  orginally  desired  to  represent  the  missile  by  an  elec¬ 
trical  analog  systen  for  the  purpose  of  organizing  and 
simplifying  the  matheuatical  analysis  Involved. 

It  would  be  ideal  if  this  analysis  could  be  performed  on  a 
purely  mathematical  basis.  This  is  possible  in  the  linear 
case.  Here  the  Laplace  Transform  is  applicable  and  matrix 
methods  exist  to  handle  the  simultaneous  equations  (2-13), 
(2-lU)  and  (2-19).  Unfortunately,  the  Laplace  Transformation 
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is  a  linear  traneornatieB  and  is  not  applicable  to  problems 
of  this  type.  Furthermore,  the  utility  of  the  Laplace 
transformation  is  destroyed  when  the  coefficients  of  the 
equations  are  functions  of  time. 

There  does  not  appear  to  be  any  uniform  mathematical  procedure 
for  handling  non-linear  cases  which  would  be  the  counterpart 
of  the  Laplace  Transform  for  linear  cases,  references  1,  8, 

9,  10.  The  state  of  the  art  in  non-linear  analysis  does  pro- 
ride  means  of  solving  certain  types  of  ordinary  non-linear 
differential  equations  in  closed  form.  In  general,  however, 
no  procedure  is  known  for  obtaining  closed  form  solutions. 

In  most  eases  it  Is  necessary  to  resort  to  iterative  procedures 
which  will  yield  approximations  to  the  solutions.  It  is  always 
possible  to  solve  a  non-linear  differential  equation  by  numeri¬ 
cal  methods.  However,  unless  the  problem  is  trivial  the  amount 
of  work  Involved  necessitates  the  use  of  a  digital  computer. 

Solutions  to  simultaneous  non-linear  differential  equations 
can  sometimes  be  obtained  by  the  process  of  elimination  of 
variables.  In  some  cases,  however,  the  nature  of  the 
equations  is  not  such  that  the  elimination  process  is 
possible.  Equations  involved  in  this  research  are  of  this 
type.  There  appears  to  be  no  method  other  then  to  resort  to 
numerical  methods  for  solving  such  equations.  It  is  common 
practice  to  analyze  systems  giving  rise  to  equations  of  this 
type  by  employing  analog  or  digital  computers. 
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The  analog  computer  le  ideally  suited  to  such  analysis  where 
extreme  accuracies  are  not  required* 

One  of  the  most  important  features  of  the  solution  of  a  non« 
linear  sysvem  is  the  stability  condition*  The  linear  system 
is  either  stable  or  unstable*  Stability  is  not  determined 
by  the  drive  function  or  initial  conditions*  In  sharp  con¬ 
trast  with  this,  however*  a  non-linear  system  may  have 
regions  of  stability  and  regions  of  Instability.  In  this 
case  stability  of  the  system  may  be  a  function  of  the  initial 
conditions  or  the  drive  ftinction*  Stability  criteria  have 
been  developed  by  Llapounoff*  reference  ^  and  Polnca]^* 
reference  10*  but  these  apply  only  to  non-linear  systems  with 
one  degree  of  freedom*  For  this  reason  the  analog  computer 
is  frequently  used  specifically  for  the  study  of  stability 
of  systems  with  multiple  degrees  of  freedom*  This  is  because 
the  analog  computer  is  Itself  a  physical  system  so  arranged 
as  to  have  the  characteristics  of  the  system  under  study* 

The* conditions  under  which  the  computer  shows  instability 
can  be  directly  translated  to  conditions  under  which  the 
system  under  study  would  similarly  becosie  unstable*, 

In  order  to  establish  the  stability  condition  of  the  missile 
system  it  is  necessary  to  obtain  the  best  possible  mathematical 
model  of  that  system*  This  implies  that  alterations  of  the 
mathematical  model  by  virtue  of  linearization  or  simplification 
or  making  approximations  can  result  in  erroneous  information 
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concerning  the  etablllty  condition*  For  this  reaaon  it  is 
considered  unwise  to  Indulge  in  such  simplification*  It  is 
therefore  evident  that  any  analog  circuit  or  simulator 
representing  the  missile  system  must  itself  contain  all 
necessary  non-linear  terms* 
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SECTIOM  III 
THE  SIMULATOR  CIRCUIT 

3.1  ChoAot  of  Compoh^ntB 

The  choice  of  componente  coaprislng  a  elnulator  circuit 
depends  prinarily  on  the  circuit  operations  to  be  performed* 
Those  involred  In  this  study  were  giren  in  Section  I  and 
are  repeated  here  for  eonrenienee: 

Amplification 
Summation 
Integration 
Inrersion  (reciprocal) 

Multiplication 

Function  Generation  (trigonometric) 

The  apparatus  capable  of  performing  each  of  these  operations 
will  be  discussed  separately* 

3*1*1  Amplification.  This  is  an  operation  on  an  input  8lg> 
nal  Inyolving  scalar  multiplication*  The  output  may  be 
posltlre  or  negatlre  and  obeys  the  eiiuation 

«o  •  i*  *1 

where  0^  is  a  generalized-  coordinate  representing  tl^e  putput 
quantity  and  0^  represents  the  input  quantity*  The  scalar 
multiplier  constant  K  may  be  greater,  equal  to,  or  less  than 
unity*  When  K  is  less  than  unity  the  operation  may  be  inter- 
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pretad  as  attanuatlon* 


A  derloe  used  to  perform  amplification  (amplifier)  Is  Intended 
to  obey  equation  (3-l)»  hoveTer'ln  practical  cases  It  does  so 
with  limitations  on  frequency,  amplitude,  gain  (K),  and  noise. 
An  Ideal  amplifier  would  obey  equation  (3-1)  without  limita¬ 
tions  and  have  infinite  Input  impedance  and ‘^ero  output  imped¬ 
ance.  Practical  amplifiers  In  general  have  the  following 
properties  t 

JL.  Moderate  Input  Impedance 
Moderate  output  impedance 
Limited  frequency  response 
d.  Limited  amplitude  range 

jo.  Variations  in  gain  with  signal  level  (nonlinearity) 

£•  Variations  in  gain  with  time 

j;.  Drift  (d.c.  amplifiers  only) 

h.  Noise  In  output 

These  properties  give  rise  to  departure  from  ideal  performance. 
The  evaluation  of  an  amplifier  depends  upon  the  specifications 
dictated  by  the  specific  application.  In  any  given  case  some 
properties  may. be  more  Important  than  others.  For  use  in  an 
analog  simulator  properties  c  and  d  offer  no  problem.  The 
important  properties  are  b,  e,  f,  g,  and  h.  Such  an  amplifier 
must  have: 

1.  Very  low  output  Impedance 

2.  Freedom  from  gain  variations  from  any  cause 

3.  Extremely  small  drift 
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lit  Lpw  noia*  output 

It  la  wall  known  that  all  of  thaaa  eharacterlatiea  can  ba 
obtalnad  by  tha  uaa  of  negatira  faadbaok  prowidad  tha  orarall 
gain  ia  low  (In  tha  ordar  of  1  to  20).  In  order  to  prorida 
iaolatlon  tha  output  inpadanea  nuat  be  rery  low.  Tha  uaa  of 
nagatira  feedback  aceonpliahaa  thla  and  at  tha  aane  tina 
randera  tha  amplifier  immune  to  iaflueneea  which  cauae 
variation 'a  in  gain.  Freedom  from  drift  la  accompllahed  by 
aubaldlary  a tablliiatlon  cireuita*  Packaged  ampllfiera  for 
thia  application  are  commercially  available  and  are  neceaaary 
to  any  analog  ‘circuit  to  ba  propoaad. 


3.1«2  Summation.  Thla  la  readily  accompllahed  by  a  paaalve 
circuit  aa  shown  in  Figure  9  where  the  input  and  output 
quantltiea  are  voltage.  The  tranafer  function  for  thia  cir* 
cult  iat 


a 


B 


(3-2) 


where  la  the  Thevenin  raalatanea  at  tha  output  terminala 
with  tha  load  realatanoe  connected.  Equation  (3-2)  ahowa 
the  aummation  operation  and  reveala  the  gain  function  to  ba 
a  constant  leaa  than  unity.  The  advantagea  of  thla  circuit 
are  Ita  almpllcity  and  reliability.  One  diaadvantage  Ilea  in 
the  fact  that  it  la  difficult  to  deaign  when  the  valuaa  of 
gain  for  varloua  inputa  are  different,  that  la,  changing  a 
resistance  on  one  input  affects  the  gain  for  all  other  inputs. 
The  load  realatanoe  also  affects  the  gain;  hence,  the  circuit 
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Fig.9.  PASSIVE  SUMMATION  CIRCUIT 
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Fig.  10.  IMPROVED  SUMMATION  CIRCUIT 
USING  HIGH  GAIN  D.C.  AMPLIFIER. 
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doea  not  exhibit  the  desirable  characteristic  of  isolation* 


Presumably  the  attenuation  resulting  from  multiple  Inputs 
in  the  ease  of  the  simple  resistire  network  can  be  made  up 
by  following  the  summation  network  with  a  low  gain  stabilised 
amplifier t  But  if  this  is  to  be  done,  it  is  more  expedient 
to  incorporate  the  amplifier  in  the  summation  circuit  as 
shown  in  Figure  10  where  summation,  amplification,  and 
isolation  are  accomplished  simultaneously  yielding  the 
further  advantage  that  the  gains  for  the  various  Inputs  are 
now  Independent  as  shown  by  the  transfer  functlont 
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(3-3) 


This  circuit  is  decidedly  superior  to  the  cascaded  summation 
circuit  and  amplifier  and  is,  therefore,  commonly  employed  in 
analog  computers* 


3*1*3  Integration*  A  passive  integrator  is  shown  in  Figure 
11  working  into  a  load  resistance  R^*  The  output  voltage  is 

taken  across  R^*  To  determine  the  transfer  function  of  the 

o 

circuit  it  is  helpful  to  obtain  the  equivalent  circuit  of 
Figure  12  by  applying  Thevenin*s  theorem  at  the  capacitor 
terminajLs*  In  this  circuit  e^  and  R^  are  the  equivalent 
Thevenin  values*  In  operational  form  the  transfer  function 
is  t 
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FIG  II  CIRCUIT  OF  PASSIVE  INTE6RATI0R 


Rt 


FIG  12  EQUIVALENT  CIRCUIT  OF  PASSIVE  NTE6RATOR 
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(3-U) 


This  equation  shows  that  the  circuit  is  not  a  true  integrator 
because  of  the  appearance  of  the  "1*  in  the  denominator*  If 
the  tine  constant  R^C  is  made  extremely  large,  equation 
reduces  to  the  approxinati-ont 


•  -  «o 


.  *1 


(3-5) 


This  reveals  that  to  obtain  good  Integration  requires  severe 
attenuation  because  of  the  large  time  constant* 


Another  possible  passive  Integrator  makes  use  of  inductance 
In  a  circuit  which  Is  the  dual  of  Figure  11*  The  output 
variable  is  then  current  rather  than  voltage.  Such  a  circuit 
offers  no  advantage  from  an  operational  standpoint  but 
possesses  the  disadvantage  that  a  satisfactory  inductor  is 
extremely  difficult  to  obtain  because  of  non-linearity  and 
hysteresis  of  the  core  material*  The  most  satisfactory 
Integrator  is  the  type  which  consists  of  an  operational 
amplifier  in  a  standard  analog  integrator  configuration* 

Even  this  is  not  a  perfect  Integrator  but  it  is  many  thousands 
of  times  superior  to  the  best  passive  Integrator*  In  this  case 
the  transfer  function  ist 


where  0  is  the  open  loop  gain  of  the  amplifier  and  may  be  of 
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the  order  of  magnitude  of  10^  or  10^*  In  comparing 
equations  (3-6)  and  (3-I4)  It  can  be  seen  that  this 
obTlatee  the  difficulty  with  the  paselre  circuit. 

Attenuation  le  not  serere  since  It  is  no  longer  necessary 
to  make  the  time  constant  large.  In  addition  the  output 
impedance  la  extremely  low  thus  proridlng  Isolation, 

3,1,U  Multiplication,  This  Is  one  of  the  most  difficult 
operations  to  perform  by  electrical  means.  Since  the 
multiplier  device  has  two  Inputs  and  one  output,  it  is 
fundamentally  a  modulator.  In  a  sense  one  input  quantity 
modulates  the  other.  No  passive  device  is  known  which  can 
perform  this  operation  with  the  possible  exception  of  a 
solid  state  semiconductor  unit  employing  the  Hall  effect. 

Here  the  output  voltage  is  the  product  of  an  input  current 
and  an  input  magnetic  field  according  to  the  equation 

-  Kj^IB  (3-7) 

where  is  a  constant  determined  by  the  physical  configuration' 
of  the  semiconductor  element,  I  is  an  input  current  through 
the  element,  and  B  is  the  magnetic  flux  density  normal  to  the 
element. 

In  practical  cases  I  is  measured  in  amperes  and  B  in 
kllogausses.  Even  at  these  levels  the  output  voltage  is 
usually  less  than  a  volt  because  of  the  small  value  of 
One  of  the  major  difficulties  encountered  in  building  a 
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Hall  multiplier  le  the  problem  of  obtaining  high  magnetic 
flux  denelty*  This  requires  the  use  of  a  magnetic  core 
exhibiting  good  linearity  and  low  hysteresis* 

In  application  the  Hall  multiplier  would  be  arranged  to 
hare  Toltage  Input  and  output*  Thus  It  would  be  necessary 
to  have  I  proportional  to  e^^  (input  No*  l)  and  B  proportional 
to  (input  No*  2)*  Because  of  the  low  output  voltage  e^ 

It  would  be  necessary  to  provide  cascade  amplifiers  to  bring 
e^  up  to  a  useable  level*  The  need  for  such  amplification 
destroys  any  advantage  of  having  a  passive  multiplier* 

Other  forms  of  solid  state  multipliers  have  been  investigated 

but  none  have  been  found  which  are  better  than  the  Hall  type* 

For  examplej  it  is  possible  to  perform  multiplication  by 

resistance  modulation.  That  Is,  a  device  whose  output  circuit 

% 

la  Isolated  from  the  Input  circuit  and  whose  output  resistance 
is  inversely  proportional  to  the  input  voltage  e^*  The  second 
Input  voltage  62  is  then  connected  across  the  output  terminals 
and  the  resulting  current  is  proportional  to  6^02*  A  practical 
unit  based  on  this  principle  requires  several  active  biasing 
circuits  which  again  destroys  the  passive  nature  of  the  device* 

A  great  deal  of  commercial  development  has  been  done  on 
function  multipliers  and  some  satisfactory  units  are 
commercially  available*  These  are  all  active  devices;  purely 
passive  multipliers  apparently  do  not  exist* 


3*1*5  Inrerslon*  This  operation  is  closely  related  to 
nultlplloatlon  because  it  constitutes  nultipllcatlon  by 
a  reciprocal.  Passive  devices  for  Inversion  are  more 
difficult  to  provide  than  those  for  multiplication*  This 
is  true  because  the 

Llm  e^ 

does  not  exist  for  an  inverter*  If  one  establishes  a  lower 
bound  on  e^,  the  operation  requires  a  finite  but  high  output 
for  small  input*  Also  the 

Llm  ■  0 

must  be  satisfied*  No  known  passive  device  is  capable  of 
performing  this  operation*  It  is  evident  that  the  resis¬ 
tance  modulator  mentioned  in  3*l*ii  might  be  applicable  to 
this  case  provided  the  output  resistance  is  made  proportional 
to  the  input  voltage  e^  and  a  constant  voltage  source  is 
connected  to  the  output  terminals.  The  output  current  would 
then  be  proportional  to  the  reciprocal  of  the  input* 

From  a  practical  viewpoint  a  resistance  modulator  inverter 
would  require  numerous  biasing  circuits  subject  to  critical 
adjustment*  The  modulator  would  consist  of  a  light  source 
focused  on  a  photoconductive  cell*  The  photoconductive  cell 
is  electrically  isolated  from  the  input  circuit  and  its 
resistance  is  the  output  quantity*  The  natural  characteristic 
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of  a  photoeonductive  cell  ahows  a  decreasing  resistance  for 
increasing  light  flux.  In  order  to  make  the  resistance 
increase  vith  light  flux  the  input  to  the  cell  must  be 
biased  to  produce  maximum  light  intensity  with  minimum  input 
▼oltage.  In  addition,  the  natural  non-linear  characteristic 
of  the  photoconductive  cell  would  require  some  form  of  phy¬ 
sical  linearizing  technique.  Such  devices  also  exhibit  a 
drift  characteristic  which  is  not  necessarily  serious  when 
performing  multiplication  but  probably  would  be  serious  when 
performing  inversion.  A  small  drift  would  repreujnt  a  large 
percent  change  in  input  level  for  small  values  of  e^. 

Inversion  can  be  successfully  accomplished  by  most  commercial 
multipliers j  in  fact  many  have  •multiply-divldo“  switches  to 
allow  either  operation. 

3.1.6  Function  Qeneration.  This  is  an  operation  involving 
the  trigonometric  functions  "sine"  and  •cosine*  where  the 
function  generator  would  be  required  to  provide  output 
voltages  representing  the  sine  and  cosine  of  the  input 
voltage  over  a  range  of  plus  or  minus  two  quadrants  (0  to 
180®  in  either  direction).  It  is  possible  to  generate  these 
functions  by  passive  resistance  diode  networks.  The  pro¬ 
cedure  for  the  design  of  such  networks  is  given  in  nearly  all 
texts  on  analog  computers.  It  will  be  found,  however,  that 
to  cover  two  quadrants  will  require  a  separate  generator  for 
each  quadrant  with  all  outputs  combined  in  a  summation  circuit. 


The  cosine  generator  will  require  biasing  one  input  of  the 
suBunauiOB  circuit  in  order  to  have  unity  output  with  aero 
input.  The  overall  function  generator  circuit  cannot  be 
coapletely  passive  because  two  of  the  input  voltages  to  the 
sunmatlon  circuit  will  have  to  be  Injected  in  the  negative 
polarity  in  order  to  provide  subtraction.  Changing  sign 
cannot  be  done  passively  and  requires  an  operational  amplifier. 
In  this  event  it  will  be  simpler  to  combine  the  operational 
amplifier  with  the  resistance-diode  network  to  produce  an 
active  function  generator  requiring  fewer  component  parts. 

Such  generators  are  very  reliable  and  once  adjusted  to  pro¬ 
duce  the  desired  function  do  not  require  frequent  readjust- 
mer.  t. 

The  need  for  function  generators  of  the  type  under  discussion 
is  such  a  common  occurrance  in  the  field  of  analog  analysis 
that  complete  packaged  units  are  commercially  available  and 
require  only  an  initial  adjustment  in  order  to  obtain  the 
desired  function.  In  the  case  a-t  hand  one  unit  would  be 
required  for  the  sine  function  and  a  separate  unit  for  the 
cosine  function.  The  accuracy  with  which  the  function  can 
be  generated  depends  upon  the  number  of  segments  used  and  can 
be  made  quite  adequate  for  use  in  the  simulator  under 
consideration. 
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3*2  Choice  of  Circuit 


3*2«1  Apparatus.  It  was  shown  In  Section  II  that,  because 
of  the  unilateral  coupling  exhibited  by  equations  (3-1) 
through  (3-3),  It  Is  not  possible  to  simulate  the  missile 
system  by  a  purely  passive  electrical  circuit.  The  unilateral 
coupling  terms  require  the  employment  of  active  Isolation 
devices.  The  electronic  amplifier  with  Inverse  feedback 
(operational  amplifier)  is  ideally  suited  to  this  application* 
The  Isolation  amplifier  must  satisfy  all  the  requirements  out¬ 
lined  in  Section  3*1*1* 

Figure  7,  which  shows  the  functional  arrangement  of  the 
simulator,  reveals  that  many  amplifiers  will  be  needed.  In 
this  case  it  would  be  desirable  to  have  all  amplifiers  alike 
and  interchangeable  so  that  one  or  two  spare  amplifiers  will 
suffice  for  maintenance.  This  in  turn  suggests  an  analog 
operational  amplifier  capable  of  operating  at  zero  frequency* 
Commercial  units  with  chopper  stabilization  are  available 
from  many  manufacturers.  Even  if  the  circuit  of  Figure  5 
were  to  be  used  to  represent  the  missile  system,  it  would  be 
necessary  to  use  six  such  operational  amplifiers. 

Now- that  it  has  been  established  that  amplifiers  are  necessary 
for  the  purpose  of  providing  isolation  it  is  clear  that  they 
can  simultaneously  be  used  for  amplification  wherever  necessary. 
It  is  now  expedient  to  consider  the  three  operations!  summation. 
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int«gration|  and  function  generation*  These  are  the 
operations  irhich  can  be  perfomed  by  passive  circuits.  How¬ 
ever,  as  shown  in  Section  3*1  these  circuits  do  not  possess 
the  property  of  isolation  and,  therefore,  must  be  cascaded 
with  anpllflers  for  that  purpose.  It  is  better  circuit  design 
and  perfomanee  is  superior  when  the  passive  circuit  is  eon- 
bined  with  the  anpllfier  as  outlined  in  any  standard  text  on 
analog  cewputlng  circuits. 

In  considering  aultipllers  and  dividers  there  appears  to  be 
no  choice  but  to  resort  to  the  use  of  comnerclal  units  which 
are  designed  specifically  for  performing  these  operations  with 
a  high  degree  of  accuracy* 

3*2.2  Theoretical  Considerations*  It  was  shown  in  Section  2 
that  because  of  the  nature  of  the  equations  involved,  it  is 
not  possible  to  establish  a  simple  straight  forward  nathematical 
procedure  for  analyzing  the  missile  system  in  terms  of  the 
analog  system.  .  This  is  primarily  because  the  state  of  the  art 
in  non-linear  analysis  does  not  provide  for  this  possibility* 
Although  the  analysis  of  linear  systems  has  been  highly 
developed  and  many  mathematical  procedures  such  as  Laplace 
Transform  are  readily  available,  the  same  Is  not  true  with 
non-llnear  systems  particularly  with  more  than  one  degree  of 
freedom*  There  does  not  even  appear  to  be  a  single  procedure 
for  the  solution  of  non-llnear  equations  representing  one 
degree  of  freedom.  The  situation  la  further  complicated  by 


the  fact  that  the  stability  problem  in  non-linear  analysis 
la  much  more  difficult  than  for  linear  systems.  Although 
stability  criteria  have  been  developed  for  special  cases 
there  is  no  general  stability  criteria  for  non-linear 
systems  as  such* 

In  the  absence  of  a  mathematical  procedure  it  then  becomes 
important  that  the  analog  circuit  allow  analysis  of  the 
system  in  terms  of  the  response  of  the  analog  system.  This 
is  a  common  and  effective  method  of  engineering  analysis. 

In  the  light  of  this  it  is  Important  that  the  analog  circuit 
be  flexible  to  permit  modifications  to  coincide  with 
possible  modifications  of  the  behavior  equations.  It  is  also 
important  that  the  circuit  be  such  that  it  is  convenient  for 
setting  coefficients  and  introducing  initial  conditions. 

One  very  Important  factor  is  that  it  must  permit  use  of 
functions  for  which  equations  may  not  be  known.  In  the  case 
at  hand  this  would  be  wind  profile  drives.  This  can  readily 
be  accomplished  in  an  analog  circuit  by  use  of  a  standard 
curve  follower  device  which  is  commercially  available. 

An  analog  circuit  arranged  to  acc\irately  represent  the  physical 
system  will  permit  a  study  of  the  stability  uondition  by 
observation  rather  than  by  theoretical  means*  The  information 
gained  by  this  procedure  is  an  extremely  Important  part  of  the 
analysis;  in  fact,  it  is  because  of  this  that  approximations 
of  the  original  equations  could  not  be  made. 


3.2.3  Recomnended  Circuit.  In  view  of  the  dlecueelon 
covered  in  Sec.  3.1.1  end  3.1.2  it  le  recommended  that  the 
circuit  of  Figure  8  be  used  for  the  analog  of  the  mleelle 
ayetem.  This  circuit  la  compriaad  of  atandard  commercial 
components  all  of  vhlch  possess  the  property  of  isolation* 

It  affords  ease  in  setting  coefficients  and  insertion  of 
initial  conditions.  It  provides  access  for  reading  out  all 
Important  variables  and  input  points  for  insertion  of  wind 
profile  drive  functions.  It  allows  easy  modification, 
should  this  become  necessary,  and  will  perform  with  satis¬ 
factory  accuracy  and  reliability. 

This  circuit  was  chosen  because  it  affords  simplicity  of 
operation  and  maintenance.  Other  circuits  were  considered 
but  none  were  simpler  or  contained  fewer  circuit  components. 


3.2, U  List  of  Components, 

Item  No,  Description  Quantity 

1  Operational  Amplifier  20 

2  Function  Multipliers  12 

3  Curve  Followers  U 

k  Function  Generators  6 


This  list  covers  all  needed  functional  devices  allowing  for 
modification  and  spares.  In  addition  to  these  items  it  will 
be  necessary  to  furnish  a  cabinet,  circuit  junction  board, 
and  other  minor  accessories  according  to  the  wishes  of  the 
group  using  the  simulator. 
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It  may  appear  that,  since  the  components  listed  are  commonly 
used  on  analog  computers.  It  would  be  possible  to  purchase  an 
analog  computer  complete  as  such.  This  may  be  possible; 
however,  the  array  of  devices  listed  Is  not  usually  found  on 
a  standard  machine.  Further  a  removable  patch  board  Is  not 
essential.  This  decision  should  be  made  by  appropriate 
personnel  at  AOHC. 
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SECTION  IV 

CONCLUSIONS  AND  RECOMMENDATIONS 
U.l  ConclttBlona 

It  l8  eoncludad  from  this  rosoareh  thatt 

a.  It  was  not  found  possible  to  represent  the  missile 
system  by  a  completely  passire  electrical  analog 
circuit* 

b«  Due  to  the  severe  non-linear  character  of  the  set 
of  behavior  equations  used  to  describe  missile 
dynamics.  It  is  not  possible  to  find  a  mathematical 
procedure  for  analysis  which  would  be  the  counter¬ 
part  of  the  Laplace  Transform  for  linear  analysis. 
£,  Because  of  the  importance  of  the  stability  problem 
in  non-linear  analysis,  it  is  essential  to  use  the 
most  accurate  mathematical  model  possible  for  the 
study  of  missile  dynamics.  Thus;  linearization 
or  other  simplifying  techniques  should  not  be 
applied  to  the  mathematical  model, 
d.  The  term  "'indlcial  admittance"  is  a  property  of  a 
linear  system  and  has  significance  insofar  as  it 
is  used  in  conjunction  with  the  superposition 
integral  or  other  linear  analytic  techniques. 

Since  superposition  does  not  apply  to  non-linear 
systems,  indlcial  admittance  as  such  does  not  have 
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equivalent  algniflcanee  for  non>llnear  syeteme 
and,  therefore,  for  the  missile  system* 

Analysis  of  missile  dynamics  can  probably  best 
be  accomplished  utilizing  an  analog  simulator  and 
interpreting  the  missile  response  in  terms  of  the 
simulator  response*  This  amounts  to  using  the  sirnu* 
lator  to‘  bypass  the  mathematical  difficulties 
arising  out  of  the  required  non-linear  analysis* 

li.2  Recommendations 

As  a  result  of  the  research  covered  by  this  report  the 
following  recommendations  are  offered: 

^*  An  analog  simulator  should  be  built  or  purchased 
as  outlined  in  Section  3*2.3* 

^*  Any  contemplated  program  for  the  development  of 

special  analog  devices  should  be  abandoned*  This 
is  because  the  product  of  such  development  does 
not  promise  to  be  superior  to  the  circuit  recommended,  . 
either  from  the  standpoint  of  circuit  performance  or 
information  it  would  reveal  to  personnel* 

£*  Further  research  should  be  carried  out  directed  to¬ 
ward  development  of  useful  mathematical  techniques 
for  handling  systems  of  non-linear  differential 
equations,  preferably  of  the  type  encountered  in 
missile  studies*  Such  research  should  be  conducted 
by  a  group  consisting  of  applied  mathematicians  and 
engineers  and  should  be  planned  as  a  long  term  program. 
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